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A solar dish Stirling cogeneration system is considered to provide energy demands of a residential build¬ 
ing. As energy demands of the building and output power of the engine are functions of weather condi¬ 
tion and solar irradiation flux, the benchmark building was considered to be located in five different cities 
in Iran with diverse climatic and solar irradiation conditions. The proposed solar dish Stirling micro-CHP 
system was analyzed based on 3E analysis. The 3E analysis evaluated primary energy saving analysis 
(energy analysis), carbon dioxide emission reduction (environmental analysis) and payback period for 
return of investment (economic analysis) and was compared to a reference building that utilized primary 
energy carriers for its demands. Three scenarios were considered for assessment and sizing the solar dish 
Stirling engine. In the first scenario, size of the solar dish Stirling engine was selected based on the lowest 
annual electric power demand while, in second, the highest annual electric power consumption was con¬ 
sidered to specify size of the engine. In the third scenario, a solar dish Stirling engine with constant out¬ 
put capacity was considered for the five locations. It was shown that implementing the solar dish Stirling 
micro-CHP system had good potential in primary energy saving and carbon dioxide emission reduction in 
all scenarios and acceptable payback period for return of the investment in some scenarios. Finally, the 
best scenario for selecting size of the engine in each city was introduced using the TOPSIS decision mak¬ 
ing method. It was demonstrated that, for dry weather, the first scenario was the best while, for hot and 
humid cities and temperate and humid places, the second and third scenarios were the best options, 
respectively. 

© 2013 Elsevier Ltd. All rights reserved. 


1. Introduction 

In the last half century, global energy consumption has been in¬ 
creased very rapidly and its growth is expected to continue over 
the next 50 years. World market energy consumption is projected 
to increase by 44 percent from 2006 to 2030 [1 ]. Solar power is one 
of the fastest growing sources of the renewable energy worldwide. 
Conversion of solar energy into electric power and/or thermal en¬ 
ergy for heating/cooling leads to significant reduction in primary 
energy consumption and greenhouse gas emission, especially in 
large scale applications. 

Solar thermal power systems use the heat produced by collec¬ 
tors which concentrate solar energy on the focal point of the collec¬ 
tor in order to heat up a working fluid in a heat engine/generator to 
produce power. Among three technologies of solar thermal power 
such as tower, trough and dish, dish Stirling systems have demon¬ 
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strated the highest efficiency [2]. For solar electric generation in 
the range of 1-100 kWe, Stirling engine is considered to be the 
cheapest one [3]. Although efficiency of Stirling engine may be 
low, reliability is high and cost is low. Moreover, simplicity and 
reliability are the keys to a cost-effective Stirling solar generator 
[3]. Kong et al. [4] presented an energy and economic feasibility 
analysis of a tri-generation unit using Stirling engine that utilized 
waste heat. They found that this CCHP unit could save more than 
33% of primary energy compared to a conventional stand-alone 
system. 

Abbas and Boumeddane [5] performed technical and economic 
analysis on a 100 MWe solar dish Stirling farm in three different 
sites. They concluded that levelized cost of the energy varied in 
the range of 11.55-23.55 C$ kW h _1 depending on the location. 
Paepe et al. [6] performed an investigation of the micro-CHP using 
three different technologies including gas engine, fuel cell and Stir¬ 
ling engine. They realized that all CHP installations, if well-sized, 
could lead to reduction in primary energy saving and carbon diox¬ 
ide emission. They also showed that gas engines had the best per¬ 
formance. Furthermore, their economic analysis demonstrated that 
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Nomenclature 


A area (m 2 ) 

C purchased cost (US $) 

c energy unit cost (US $ kW h -1 ) 

CDE carbon dioxide emissions (kg h -1 ) or (kg year -1 ) 

CDER carbon dioxide emissions reduction 

DM direct normal irradiation 

E electric power (kW) 

E b electric power used by hot water boiler (kW) 

E ch electric power used by absorption chiller (kW) 

/ supplied fuel energy (kW h) 

HVAC heating, ventilating and air conditioning 
H.W.B hot water boiler 

NG natural gas 

Q. thermal energy (kW) 

Qcc thermal energy delivered to cooling coils (kW) 

Q cb thermal energy required for generator of absorption 
chiller (kW) 

Q bc thermal energy delivered to heating coils (kW) 

Q LSb thermal energy lost from hot water boiler (kW) 
QLLsv.w.m thermal energy lost from D.W.H 1 (kW) 

QLls,d.w.h2 thermal energy lost from D.W.H 2 (kW) 

QLls,h.e\ thermal energy lost from H.E 1 (kW) 

Qls,s.e thermal energy rejected from the cold-end of Stirling 
engine (kW) 

Q r thermal energy recovered at the cold end of Stirling en¬ 

gine (kW) 


T temperature (K) 

Greek symbols 

r] efficiency, useful output divided by input amount 

Subscripts 
b boiler 

BL book life 

c cooling 

cc cooling coil 

ch absorption chiller 

COLEC collector 

grid electricity supplied by grid 

he heating coil 

hw hot water 

hws domestic hot water 

in sun irradiation 

H.E heat exchanger 

L levelized 

LS heat loss 

r heat recovery 

rec heat received from the sun 

req electric energy demand 

SE Stirling engine 

sold electricity sold to grid 


fuel cells were too expensive in comparison to gas engines which 
had small rate of return for capital investment. From the economic 
aspect, installing a CHP in a residential building is not interesting 
because of the high capital cost and long period for return of the 
capital investment. However, it is required to mention that their 
Stirling engine is an externally gas burned system. Among four 
main primary movers of the micro-CHP system including internal 
combustion engine, micro-gas turbine, fuel cell and Stirling engine, 
there is only Stirling engine that can be integrated to the solar en¬ 
ergy. Therefore, it is expected to have more primary energy saving 
and more significant reduction in greenhouse gas emission in the 
case of CHP system with solar dish Stirling primary mover. Nosek 

[7] introduced a new lOkWe hybrid solar dish Stirling which 
was used in micro-CHP applications. They applied a natural gas 
burner to heat up the Stirling engine’s hot-end when electric out¬ 
put of the system was lower than 10 kWe during 24 h of the sys¬ 
tem operation in the hybrid mode. Results showed that the 
annual amount of primary energy consumption decreased to 
49 MW h. Also, the total amount of annual energy saving in the so¬ 
lar mode was 50.38 MW h, which was equivalent to 4800 m 3 of the 
natural gas consumption in the conventional gas burner. Li et al. 

[8] performed a real test on a Stirling engine driven by exhaust 
waste gas. Their test results showed that the output shaft power 
could reach 3476 W at shaft speed of 1248 RPM. It means that, 
the Stirling engines driven by mid-high temperature waste gases 
were capable of achieving valuable power for engineering applica¬ 
tions. Bravo et al. [9] made a comparison in terms of environmental 
effect between a photovoltaic power generation system and a dish- 
Stirling technology based on C0 2 -equivalent emissions. The results 
demonstrated that level of environmental impacts was similar for 
both technologies. Magri et al. [10] investigated an energetic and 
environmental analysis on a typical micro-CHP driven by a natural 
gas fired Stirling engine. Then, operation of the mentioned plant 
was compared with that of the condensing boiler to meet heating 
demands. The result showed that the micro-CHP unit can provide 


reduction in the primary energy demand with respect to a conven¬ 
tional boiler (-27.3%) and also with respect to a condensing gas 
boiler (-9.1%). Furthermore, they concluded that utilizing micro- 
CHP unit may lead to remarkable reduction of C0 2 emission. Har- 
rod et al. [11] performed an analysis on biomass-fired Stirling en¬ 
gine in order to determine how an individual component of the 
micro-CHP system such as prime mover, heat recovery system, 
auxiliary boiler, absorption chiller and heating unit coil affect the 
system’s primary energy consumption (PEC) and operational cost. 
They recommended engine size of 8 kW as the best choice so as 
to achieve the best cost of system applied for the 511m 2 area facil¬ 
ity located in Atlanta. There are a lot of valuable investigations in 
line with 3E analysis. Mago and Chamra [12] performed an analysis 
and optimization on a typical CCHP system considering energy and 
economic and environmental aspects. They introduced three crite¬ 
ria which were primary energy saving (PE), emission reduction 
(ER) and operational cost (OC) for optimization of their plant. Also, 
they applied three operational strategies following electric load 
(FEL), following thermal load (FTL) and hybrid electric thermal load 
(HETS) for performance evaluation of their proposed system. In or¬ 
der to obtain a better view of optimized value, the authors com¬ 
pared their CCHP system with a reference building which used a 
conventional HVAC system for meeting energy demands of the 
simulated building. The results explicated that, for the evaluated 
cities, the best optimization criterion was PE-0 (primary energy 
optimization) because of minimizing value of PEC while decreasing 
values of cost and CDE. Cho et al. [13] implemented an assessment 
on CCHP systems based on operational cost, primary energy con¬ 
sumption and carbon dioxide emission using different cities with 
various climatic conditions. Their results showed that CCHP sys¬ 
tems augmented the site energy consumption (SEC) with respect 
to the reference building. Results for various site locations clarified 
the fact that there was not a common trend among three optimiza¬ 
tion criteria which were primary energy optimization (PE-O), cost 
optimization (Cost-O) and carbon dioxide emission optimization 
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(CDE-O). This can be attributed to the variation of building energy 
loads, electricity and fuel costs and C0 2 emission factor for the 
evaluated cities. Alanne et al. [14] performed an evaluation and 
optimization of Stirling engine in a micro-CHP system for a resi¬ 
dential building. They searched an optimally operated configura¬ 
tion for integration of a Stirling engine-based micro-CHP system 
in a residential building by comparing performance of various sys¬ 
tems and strategies of a reference system. The results demonstrate 
that an optimally operated micro-CHP including heat recovery and 
suitable thermal storage would result in 3-5% reduction in primary 
energy consumption and C0 2 emissions when compared to a con¬ 
ventional hydronic heating system. Furthermore, their suggested 
system ensured creating annual saving in all combinations of elec¬ 
tricity and fuel price of between 0.05 and 0.15 € kW h -1 . Dorer and 
Weber [15] carried out an energy and C0 2 assessment of residen¬ 
tial micro-CHP system using a transient model in TRNSYS simula¬ 
tion program. In order to find the best cogeneration technology in 
terms of primary energy saving and C0 2 emission, various CHP 
technologies such as natural gas fueled solid oxide (SOFC), polymer 
electrolyte membrane fuel cells, Stirling and internal combustion 
engines were taken into account. The results were compared with 
a conventional building which used a natural gas-fired boiler and 
electricity from grid to meet energy demands of the building. Also, 
a ground-coupled heat pump system was analyzed for further 
comparison. The analysis demonstrated that the maximum reduc¬ 
tion of non-renewable primary energy with the evaluated cogene¬ 
ration system was 14%. In the case of C0 2 -eq emissions, most 
cogeneration systems offered reduction in the case of using Euro¬ 
pean electricity mix (up to 22%). For the case of electricity genera¬ 
tion by combined cycle power plant, maximum C0 2 -eq emission 
reduction was resulted in heat pump systems (up to 29%). More¬ 
over, the maximum decrease of C0 2 -eq emission was obtained in 
the internal combustion engine when being applied in a single¬ 
family house. E. Cardona et al. [16] performed 3E analysis and opti¬ 
mization on internal combustion engine-driven CCHP integrated to 
a reversible heat pump. The obtained results indicated that a com¬ 
plete matching between “economic” and “energy-environmental” 
choices was not feasible because of requiring “ad hoc” measures 
for each case study. 

In this paper, a typical solar dish Stirling unit was integrated by 
an HVAC system to meet electric, heating, cooling and hot water 
demands of a typical 5-storey residential building. A benchmark 
building was assumed to be located in five different cities with di¬ 
verse climatic conditions (including temperate and dry, cold and 
dry, hot and humid, temperate and humid and hot and humid 
weathers) to obtain different electric, heating, cooling and hot 
water load profiles. The proposed CHP system was evaluated based 
on 3E analysis. The 3E analysis considers energy, environmental 
and economic analysis simultaneously. An energy analysis was 
performed to illustrate benefits of the proposed system. With cal¬ 
culating site-source energy conversion factor for electricity and 
natural gas, the primary energy saving (PES) was obtained for the 
plant. Moreover, by calculating a factor that represented carbon 
dioxide emission (namely, carbon dioxide emission factor), ability 
of the proposed CHP system for reducing the carbon dioxide emis¬ 
sion was assessed in the aforementioned five different climatic 
conditions. Furthermore, the economic analysis based on the total 
revenue requirement (TRR) method was performed to evaluate 
payback period for return of investment used for providing CHP 
system. Energy, environmental and economic analysis of the pro¬ 
posed micro-CHP system with a dish Stirling prime mover was per¬ 
formed in three scenarios. In the first scenario, size of dish Stirling 
engine in each city was selected based on the minimum annual 
electric load. In this scenario, difference between the actual elec¬ 
tricity demand and output electric power of the solar dish Stirling 
engine was assumed to be delivered by the grid electricity. In the 


second scenario, the prime mover was sized based on the highest 
annual electric power demand. In the third scenario, a Stirling en¬ 
gine with constant output power capacity of 10.0 kW at maximum 
solar irradiation (1000 W m -2 ) was utilized for the proposed CHP 
system in each city. Finally, the best option for sizing the prime 
mover of the CHP system in each location was suggested using 
TOPSIS decision making method. 

2. System description 

A micro-CHP system with solar dish Stirling engine as prime 
mover as indicated in Fig. la was considered for providing electric, 
heating and cooling demands of a residential building. This micro- 
CHP system was considered to support electric, heating and cool¬ 
ing loads of a benchmark 5-story 10-flat residential building. The 
building was aligned with north with the total living area of 
960 m 2 . Complimentary information regarding the proposed build¬ 
ing is given in Fig. lb and c and Table 1. 

A hot water natural gas boiler (H.W.B) and a single effect 
absorption chiller (S.E.A.C) were utilized to supply heating and 
cooling demands of the building. Moreover, a typical free-piston 
Stirling engine integrated to a solar concentrator was provided to 
supply a part of electricity demand. The rest of electricity demand 
was provided by the local grid in the case that electricity demand 
of the building was larger than the output power capacity of the 
engine. A jacket heat exchanger (H.E 2) recovered rejected heat 
at the cold-end of Stirling engine. This recovered heat was ab¬ 
sorbed by the engine’s cooling water in a heat exchanger, namely 
D.W.H 1, in Fig. la; hence, it preheated the inlet cold tap water 
(from city water network) used for providing the hot water de¬ 
mand. The preheated water came from the D.W.H 1 flows through 
another heat exchanger (D.W.H 2), which was warmed up by the 
hot water boiler. If output power of dish Stirling system was higher 
than the electricity demand, the excessive electricity was sold to 
electricity grid. Furthermore, as is observed in Fig. la, hot water 
boiler provided heating demand of the building through heat ex¬ 
changer H.E 1 and heating coils located inside the building (nomi¬ 
nated as H.C in Fig. la). Moreover, H.W.B provided heating load of 
the generator of a single effect absorption chiller (S.E.A.C). The 
S.E.A.C provided cooling demands of the building through the room 
cooling coils (C.C). 

3. Modeling 

3.1. Estimating energy demand 

Energy demand of the building highly depended on climatic 
conditions. In order to evaluate effectiveness of the proposed mi¬ 
cro-CHP system with the solar dish Stirling engine as prime mover 
for reducing primary energy consumption and carbon dioxide 
emission when utilized in various climatic conditions, five cities 
in Iran with diverse climatic conditions were selected to obtain an¬ 
nual profiles of buildings’ electric, heating and cooling loads. These 
cities were Tehran (selected as the representative of temperate and 
dry weather), Tabriz (selected as the representative of cold and dry 
weather), Bandar Abbas (selected as the representative of hot and 
humid weather), Rasht (selected as the representative of temperate 
and humid weather) and Yazd (selected as the representative of 
hot and dry weather). Climatic outdoor design conditions of these 
cities are indicated in Table 2. 

To calculate buildings’ heating and cooling loads, the proposed 
benchmark residential building was simulated in Carrier-HAP3.4 
software (heating load also included domestic hot water load). 
Electricity demand of the reference building was classified to two 
categories including variable and fix loads. Variable loads denoted 
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Fig. 1. (a) Schematic of the solar dish Stirling micro-CHP system; (b) West view of the proposed building; and (c) Plan of the building flats. 


refrigerator loads, lightening and HVAC plant equipment electric 
consumption, varying hourly in the whole month of a year for 
the selected cities. All the electrical loads, except variable loads, 


considered the fixed loads. It is momentous to mention that fix 
loads refer to electric loads that do not vary through the seasons 
of the year. As a matter of fact, they can change from 1 h to another 
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Fig. 1 ( continued ) 


during a day. All appliance loads, except refrigerator load, were 
deemed as fix loads. Electric loads consumed by HVAX equipment 
were obtained from vendors catalogues of these equipment. Table 3 
shows the electric appliances assumed to be used in each flat of the 
building. Figs. 2-4 demonstrate heating, cooling and electricity 
load profiles of the proposed benchmark building, respectively. In 
these figures, the magnitudes assigned to each month were an 
hourly average through 24 h per day multiplied by the number 
of days in a month. 

3.2. Energy analysis of solar dish Stirling engine 

According to the BSR company record in EPRI report [17], real 
efficiency of their typical solar powered Stirling engine was re¬ 


ported as 50% of the Carnot efficiency. The Carnot efficiency is cal¬ 
culated as: 

dCarnot = ^ ~ (Jc/T'h) 0 ) 

where rjcamot is Carnot efficiency of the Stirling engine, T c is temper¬ 
ature at the cold-end of the Stirling engine (K) and T h is temperature 
at the hot-end of the Stirling engine (K). T c is assumed equivalent to 
the ambient temperature. In order to estimate hot-end temperature 
in different solar irradiation, T h = 800 °C was considered at maxi¬ 
mum DNI = 1000 W m" 2 [18] and, with a simple proportion with 
various DNIs for the selected cities (Fig. 5), values of T h were deter¬ 
mined as indicated in Table 4. Also, values of the ambient temper¬ 
ature (T a ) are given in Table 5. Values of the estimated real 
efficiency of the Stirling engine in the proposed cities during a year 
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Table 1 

Description of the simulated building in Carrier, HAP 3.4. 


Orientation 

Aligned with north 

Building type 

Residential apartment 

Number of floor 

5 

Number of flats in each floor 

2 

Area of each flat 

96.0 m 2 (7.1 m x 13.5 m) 

Average ceiling height 

3.0 m 

External wall overall (J-value 

2.179 (W/m 2 K) 

Internal wall overall (J-value 

2.839 (W/m 2 K) 

Windows overall (J-value 

3.591 (W/m 2 K) 

Roof overall (J-value 

0.684 (W/m 2 K) 

Floor overall (J-value 

0.454 (W/m 2 K) 

Total occupants 

40 

Total glass area 

120 m 2 


Table 3 

List of the considered electric appliances in each flat of the building. 


Electric 

appliance 

Number in each flat of the 
building 

Consumed power for each 
unit (W) 

Iron 

1 

1000 

LED television 

1 

150 

Hair dryer 

1 

1000 

Ventilation fan 

2 

30 

Washing 

1 

500 

machine 

Laptop 

1 

80 

computer 

Freezer / 

1 

350 

refrigerator 

Lightening 

8 

30 

lamp 


are presented in Table 6. Using Iran’s meteorological data [19], the 
number of hours during which there was solar irradiation during 
each month of the year was obtained and used to calculate primary 
energy saving and carbon dioxide emissions reduction, represented 
in Sections 2.5 and 2.6, respectively. These data are presented in 
Table 7. 

The power generated by the solar dish Stirling engine can be 
estimated based on the following expression: 

Ese = RF ■ HLF ■ A ■ f] SE ■ Qsolarjrr (2) 

where t] S e is thermal efficiency of the Stirling engine as indicated in 
Table 6 and Qsoiarjn is solar direct normal irradiation (W m -2 ) as 
indicated in Fig. 6. Moreover, in Eq. (2), Ese , RF and HLF are real out¬ 
put electric energy of the Stirling engine, reflectance factor for the 
solar dish Stirling system [7] and heat loss factor of the solar dish 
Stirling system [7], respectively. Values for RF and HLF were consid¬ 
ered 0.90 and 0.74, respectively. A in Eq. (1) is the area of the solar 
dish. The output power as the nominal size of the solar dish was 
specified in three scenarios in this research: 

• Scenario I: In this scenario, solar dish Stirling engine was sized 
based on the base electricity demand (minimum annual electric 
power demand) in the electric load profile (Fig. 4). Therefore, in 
this scenario for each location, a specific size of engine as indi¬ 
cated in Table 8 was obtained. 

• Scenario II: In this scenario, solar dish Stirling engine was sized 
based on the maximum annual electric power demand (see 
Fig. 4). In this scenario for each location, a specific size of engine 
was obtained, as indicated in Table 8. 

• Scenario III: An engine with fixed output power capacity was 
used in all locations (the proposed five cities). In this scenario, 
a solar dish Stirling engine with 10 kW nominal power obtained 
at maximum solar irradiation of 1000 W m -2 was used as prime 
movers of the micro-CHP system in all the studied sites. 

Moreover, the associated dish diameters for the considered 
sizes of Stirling engine in the selected cities are presented in 


Table 8. If output capacity of the dish Stirling system was not 
sometimes adequate to cover the whole electric demand, the dif¬ 
ference between actual electricity demand and output electric 
power of the solar dish Stirling engine would be provided by 
the electricity local grid. Therefore, the following expression can 
be given: 



0 when E se > E m 

Ereq Ese when Ese ^ E req 



where E req is electric power demand of the building (kW) and E grid is 
purchased electric power from the local grid (kW). 

Furthermore, in the case that electricity demand was less than 
the output capacity of the engine, the excess electricity was sold 
to the local grid. 



0 when E se ^ E req 

Ese E req when Ese ^ E req 



where F so/d denotes excess electricity sold to the local grid. 


3.3. Energy analysis of the entire micro-CHP system 

The total first law thermodynamic efficiency for the entire mi¬ 
cro-CHP system was obtained from the following expression: 

h overallplant (Q-c "T Qh "T Qhw H - E req )/ (Qrec "T f b "T E gr [ d E so [ d ) (4) 

where rj overall,plant is overall efficiency of the micro-CHP solar dish 
Stirling system. Q. c , Qh, Qhw and E req are the building’s cooling, heat¬ 
ing, domestic hot water heating and electrical loads in kW, respec¬ 
tively. Moreover, Q r ecjb,E grid and F so/d are the supplied solar heat to 
the hot-end of the Stirling engine (obtained from Eq. (5)) in kW, rate 
of fuel energy supplied to the hot water boiler (kW), purchased 
electric power from the local grid (kW) and electricity sold to the 
local grid (kW), respectively. 


Table 2 

Climatic data of the selected cities. 


City 

Weather 

category 

Latitude 

(°N) 

Longitude 

(°E) 

Average annual relative 
humidity (%) 

Summer design 
dry bulb (°C) 

Winter design dry 
bulb (°C) 

Summer design 
wet bulb (°C) 

Winter design wet 
bulb (°C) 

Tehran 

Temperate 
and dry 

35.7 

51.4 

40.1 

38.9 

-6.7 

23.9 

-8.8 

Tabriz 

Cold and dry 

38.1 

46.3 

53.7 

33.9 

-10.8 

18.0 

-11.8 

Bandar 

Abbas 

Hot and 
humid 

27.2 

56.4 

65 

40.6 

7.5 

31.9 

5.2 

Rasht 

Temperate 
and humid 

37.2 

50.2 

71.3 

31.9 

-2.2 

25.7 

-2.9 

Yazd 

Hot and dry 

31.9 

54.4 

31.4 

40.0 

-5.3 

18.3 

-7.3 











354 


R.S. Moghadam et al./Energy Conversion and Management 75 (2013) 348-365 



Fig. 2. Monthly heating loads (Q./,) for the selected cities. 
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Fig. 3. Monthly cooling loads (Q c ) for the selected cities. 




10 1 
9- 


^ Tehran §§ Tabriz 


Bandar Abbas 


Rasht Yazd 


JAN FEB MAR APR MAY JUN JUL AUG SEP OCT NOV DEC 


Fig. 4. Monthly electrical loads ( E req ) for the selected cities. 


Terms in numerator of Eq. (4) are load data estimated and pre¬ 
sented as loads’ profiles in Figs. 2-4, respectively. Q rec is obtained 
from the energy balance equation at the solar dish as follows, 

Crec = Qsunjrr ■ A ■ RF ■ HLF (5) 

For/i, in Eq. (4): 



where Qb and r\ b are boiler heating capacity (kW) and thermal effi¬ 
ciency, respectively. Boiler efficiency along with specifications of 
some other HVAC equipment is presented in Table 9. 

For estimation of the required capacity of the hot water boiler 
(Qb), a comprehensive first law thermodynamic analysis should 
be performed for the entire micro-CHP system. Complete list of en¬ 
ergy balance equations for the entire micro-CHP plant is given in 
Table 10. 
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Fig. 5. Monthly hot water heating loads ((Q./ 1W ) for the selected cities. 
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Table 4 

The hot-end temperature (T h ) of the solar Stirling engine installed in the selected cities. 


City T h (K) 



January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

December 

Tehran 

656.5 

763.9 

857.7 

913.3 

936.2 

939.7 

930.1 

902.1 

835.8 

753.5 

652.9 

606.6 

Tabriz 

629.3 

739.9 

840.4 

904.9 

932.6 

936.6 

924.7 

893.2 

820.8 

729.7 

629.5 

581.4 

Bandar Abbas 

759.0 

848.3 

924.5 

956.4 

960.8 

955.3 

949.3 

940.9 

906.1 

842.7 

759.9 

721.7 

Rasht 

636.3 

746.4 

844.7 

906.3 

932.9 

937.1 

926.1 

894.9 

822.5 

732.8 

629.1 

586.8 

Yazd 

702.6 

802.9 

886.7 

934.7 

949.7 

947.4 

938.0 

920.0 

870.2 

795.8 

703.3 

659.3 


Table 5 

Values of the ambient temperature (T a ) for the selected cities. 


City T a (K) 



January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

December 

Tehran 

276.3 

278 

283.3 

289.8 

294.7 

300.0 

303.3 

302.2 

298.3 

291.4 

284 

278.6 

Tabriz 

271.4 

273.6 

278.3 

284.3 

289.3 

294.5 

298.9 

298.8 

294.3 

286.9 

279.7 

274.6 

Bandar Abbas 

291.1 

292 

295.5 

299.5 

303.5 

306.2 

307.2 

307.0 

305.2 

302.4 

297.2 

292.6 

Rasht 

279.6 

278.9 

281.8 

287.1 

291.7 

295.5 

298 

297.7 

295.0 

290.0 

285.4 

281.9 

Yazd 

278.8 

281 

286.3 

292.2 

297.4 

302.5 

304.3 

302.3 

298.1 

292.4 

284.9 

280.5 


Table 6 

Actual efficiency of the Stirling engine installed in the selected cities. 


City rj rea i 



January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

December 

Tehran 

0.29 

0.32 

0.33 

0.34 

0.34 

0.34 

0.34 

0.33 

0.32 

0.31 

0.28 

0.27 

Tabriz 

0.28 

0.32 

0.33 

0.34 

0.34 

0.34 

0.34 

0.33 

0.32 

0.30 

0.28 

0.26 

Bandar Abbas 

0.31 

0.33 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

0.33 

0.32 

0.30 

0.30 

Rasht 

0.28 

0.31 

0.33 

0.34 

0.34 

0.34 

0.34 

0.33 

0.32 

0.30 

0.27 

0.26 

Yazd 

0.30 

0.33 

0.34 

0.34 

0.34 

0.34 

0.34 

0.34 

0.33 

0.32 

0.30 

0.29 


Table 7 

Monthly number of hours with solar irradiation. 


City 

January 

February 

March 

April 

May 

June 

July 

August 

September 

October 

November 

December 

Annual 

Tehran 

173.6 

182.8 

205.5 

222.8 

289.5 

345.1 

349.5 

340.3 

304.0 

256.1 

194.9 

166.1 

3030.2 

Tabriz 

123.5 

146.1 

179.5 

201.3 

269.5 

336.7 

353.6 

339.5 

303.3 

231.8 

177.3 

132.2 

2794.3 

Bandar Abbas 

221.9 

210.7 

225.4 

254.9 

312.8 

307.6 

271.8 

271.2 

265 

281.9 

251.8 

226.2 

3101.2 

Rasht 

91.7 

89.2 

89.9 

123.6 

174.0 

209.7 

220.5 

180.3 

137.7 

118.7 

96.3 

89.1 

1620.7 

Yazd 

190.3 

210.6 

218.9 

242.9 

300.7 

343.5 

345.7 

345.0 

316.0 

286.4 

224.2 

199.0 

3223.2 


In the equations given in Table 10, rh, h , E and Q are mass flow 
rate (kgs -1 ), specific enthalpy (kj kg" 1 K _1 ), electric power con¬ 
sumption by the HVAC equipment (kW) and thermal capacity of 


the HVAC equipment (kW), respectively. Furthermore, in Table 10, 
t) and COP denote thermal efficiency and coefficient of perfor¬ 
mance for heating and cooling equipment, respectively. 
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Table 8 Table 9 

General specifications of the dish Stirling engine installed in the selected cities. Input data of the simulation. 


City Electrical output capacity of the Dish diameter (m) 
solar dish Stirling engine (kW) a 



Scenario 

I 

Scenario 

II 

Scenario 

III 

Scenario 

I 

Scenario 

II 

Scenario 

III 

Tehran 

4.00 

5.84 

10.00 

5.05 

6.10 

8.00 

Tabriz 

3.58 

4.70 

10.00 

4.72 

5.50 

8.00 

Bandar 

Abbas 

9.60 

9.05 

10.00 

7.50 

7.34 

8.00 

Rasht 

7.97 

5.67 

10.00 

6.98 

6.02 

8.00 

Yazd 

3.80 

5.36 

10.00 

4.88 

5.86 

8.00 


a Electric output capacity of the solar dish Stirling engine is presented based on 
its output at maximum solar irradiation flux of 1000 W m -2 for all the scenarios. 


3.4. Evaluation of the primary energy saving (PES) 

In order to estimate primary energy saving (PES) when the solar 
dish Stirling technology was used as a prime mover of the residen¬ 
tial micro-CHP system, the primary energy consumption (PEC) of 
the proposed plant was compared with a residential building 
which utilized conventional HVAC system consuming electric and 
natural gas energy in order to follow up the same loads’ profiles 
for that building. In this regard, a factor called source-site ratio 
(SSR) was utilized as introduced in the following sections. 

3.4.1. Calculation of electricity and natural gas source-site ratios 

In order to assess the primary energy consumption, it is essen¬ 
tial to have source-site energy conversion ratios. These ratios were 
defined and taken into account for the electricity and natural gas 
consumption in this study. 

According to EPA report [20], the electricity source-site ratio 
(ESSR) is calculated as the total primary energy involved in the 
electricity generation divided by the net generation minus trans¬ 
mission and distribution (T&D) losses. This calculation is summa¬ 
rized in Table 11. The annual values for the primary energy 
consumption for the power generation, net generation and T&D 
losses were extracted from the Iran Energy Bulletin and Iran energy 
flowcharts [16,17]. 

The source-site ratio for the natural gas (NGSSR) was obtained 
first by computing sum of natural gas delivered to consumers, 
pipeline, distribution use and plant use and then by dividing the 
result by the total delivery to consumers [20]. The calculated val¬ 
ues for the natural gas source-site ratios are summarized in 
Table 12. 


Description Symbol Mean value 


Boiler efficiency 

db 

0.90 

Absorption chiller’s coefficient of performance 

COP ch 

0.70 

Heat exchangers efficiency 

dH.E 

0.75 

Cooling coil efficiency 

dec 

0.80 

Heating coil efficiency 

dhc 

0.80 


Table 10 

Energy balance equations for the entire micro-CHP plant with solar dish Stirling 
engine. 

Equation number 

Energy balance equation 

1 

m 7 • (. h u - h 7 ) - Qcc = 0 

2 

Q.cc — COP ch • Qcii = 0 

3 

Q-C — dec ■ Qcc = 0 

4 

rh u ■ (hu - hu) - Qhc = 0 

5 

Qh — dhc Qhc ~ 0 

6 

Qb -fbdb = 0 

7 

m 5 ■ (h 5 -h u )-Qch =0 

8 

m 2 o • (h 2 o - h 22 ) + rh 6 • (h 6 - h 12 ) - Qch - Ech = 0 

9 

rhu • (/in - hi 5 ) + hii 6 • (hi 6 - hio) + Qls,h.ei = 0 

10 

th 2 6 • (h 2 6 - h 2 5 ) + rii28 • (h 2 8 - h 2 7 ) + Chs,Dwm = 0 

11 

mg ■ (hg - h 8 ) + rh 2 g • (h 2 9 - h 2 s) + Qls,dwh2 = 0 

12 

Q b - rh 3 ■ (h 3 - hig) = 0 

13 

Qb + m 23 ■ h 23 - m x ■ h x - m 2 ■ h 2 + Chs,b - E b = 0 

14 

QlSS.E ■ d\\.E 2 - hi 24 ■ (h 2 4 ~ tl 2 6 ) = 0 

15 

Qhw ~ tn 2 7 • (/130 - h 2 7) = 0 

16 

Qr - tn 21 ■ (h 2 s — ft 2 7) = 0 

17 

Qrec = Qsun.irr ' A ■ RF ■ HLF 

18 

Qhw ~ Qr — bl 2 8 • (/l30 - h 2 8 ) = 0 


Heat losses from surface of heat exchanger were assumed to be 1% of exchanger’s 
heat duty. For the boiler, it was assumed to be 2% of the boiler capacity. 


3.4.2. Analysis of the plant based on primary energy consumption 
The primary energy consumption for the micro-CHP solar dish 
Stirling system ( PEC CH p-soiar ) was: 

PECcHP-Solcir = ( Keq ~ ^Se) ' ESSR 

+ (fe)- <i ' +C! * + (^)) NCSSS 0> 

where rj hws and rj HE{ are hot water system efficiency and efficiency 
of heat exchanger #1 (indicated in Table 9), respectively. ESSR and 
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Fig. 6. Monthly solar irradiation ( Q. so iarm ) f° r the selected cities. 
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Table 11 

Electricity source- site ratios. 


Year 

PEC a for power 
generation (GW h) 
[21,22] 

Net generation 
(GWh) [21,22] 

T8jD b loss 
(GWh) 
[21,22] 

SSR C 

2003 

416168.03 

142942.87 

26575.41 

3.570 

2004 

460129.44 

159739.57 

30341.77 

3.550 

2005 

499531.89 

170252.98 

33086.72 

3.640 

2006 

N.A 

184589.16 

32747.60 

N.A 

2007 

523629.26 

196030.74 

45114.98 

3.470 

2008 

592502.13 

206163.71 

45518.69 

3.680 

2009 

610282.08 

212960.05 

43005.63 

3.590 

Average(2003-2009) 



3.580 


a Primary energy consumption. 
b Ransmission and distribution. 
c Source-site ratio. 


Table 12 

Natural gas source-site ratio. 


Year 

Some of pipeline, distribution use, 
plant fuel and delivery to consumers 
[20] (MBOE a ) 

Delivery to 
consumers [20] 
(MBOE) 

Source- 

site 

ratio 

2003 

546.20 

334.00 

1.635 

2004 

615.10 

387.10 

1.612 

2005 

677.60 

427.00 

1.580 

2006 

N.A 

N.A 

N.A 

2007 

777.80 

471.80 

1.648 

2008 

820.20 

476.50 

1.720 

2009 

859.90 

519.70 

1.650 

Average(2003-2009) 


1.640 


a MOBE: Million barrel oil equivalent. 


NGSSR denote electricity and natural gas source-site ratios, respec¬ 
tively. It is required to mention that, in the case that the electricity 
demand was less than the output capacity of the engine, the exces¬ 
sive electricity was sold to the local grid. Therefore, in such a case, 
(. E req - E S e ) was negative and caused reduction from the primary en¬ 
ergy consumption ( PEC CH p-soiar )• 

Furthermore, the primary energy consumption for the reference 
building with conventional HVAC system PEC re f was: 


PEC ref = ( Ereq ) ' ESSR + 


/ Qhws 

\fl hws 


+ Qch + 



.NGSSR 



Then, the primary energy saving (PES) in the case of using micro- 
CHP solar Stirling system was: 

PES = PECref - PECcHP-Solar (9) 


3.5. Evaluation of carbon dioxide emission (CDE) 

Electricity and natural gas emission factors were required for 
evaluating carbon dioxide emission for the proposed system. These 
values could be obtained from [23] as reflected in Table 13. For 
simplification purposes, it should be mentioned that these factors 
were considered the same in all the evaluated locations. 

In order to show the environmental protection effect of the pro¬ 
posed system, it was compared with the reference building, which 
utilized a similar conventional HVAC system. To this end, carbon 
dioxide emissions from the CHP-solar Stirling system ( CDE C HP-soiar ) 
were formulated as follows: 

CDEcHP-Solar = (Ereq ~ Ese) ' CDEFe 


Table 13 

Conversion factor of carbon dioxide emissions for the electricity and natural gas and 
energy price [23,31,32]. 


Energy 

carrier 

Conversion factor 
carbon dioxide 

emissions 

Purchased cost of 
the energy 

Sold price of the 
energy 

Electricity 

0.598 kg kW h -1 

0.0150US $ kW h -1 

0.4200US $ kW h -1 

Natural 

0.185 kg kWh” 1 

o.o^eusskwh- 1 

N.A 

gas 





where CDEF e and CDEF NG denote carbon dioxide emission factors for 
the electricity and natural gas, respectively. Also, the carbon dioxide 
emission for the reference building can be obtained as follows: 


CDEref = (Ereq) ' CDEF E + 


| Qhws 
y d hws 


+ Qch + 



.CDEF ng (11) 


Therefore, value of the carbon dioxide emission reduction (CDER) 
can be determined as: 


CDER — CDEref — CDEcHP-Solar (12) 

The CDER parameter showed the amount of carbon dioxide emis¬ 
sion which was diminished on account of using the micro-CHP solar 
Stirling system in comparison to the reference conventional resi¬ 
dential building. 


3.6. Economic analysis 

In this section, associated cost of the system capital investment 
for the solar dish Stirling system was compared with saving in the 
cost of primary energy consumption. In this regard, difference be¬ 
tween capital investments of the micro-CHP system and reference 
building was divided by the saving in the cost of primary energy 
saving in order to estimate payback period for each location (five 
cities in this study).The economic analysis of this research was per¬ 
formed based on the total revenue requirement method. This 
method was developed based on the procedures adopted from 
the report of Electric Power Research Institute (EPRI; 1993) [24], 
refined and presented by Bejan et al. [25]. 

The annual total revenue requirement (TRR, total product cost) 
for a system is the revenue that must be collected in a given year 
by selling all products to compensate for the system operating 
company’s all expenditure incurred in the same year and to ensure 
sound economic system operation [25]. It consists of two parts: 
carrying charges and expenses. Carrying charges are a general des¬ 
ignation for charges that are related to the capital investment 
whereas expenses are used to define costs associated with the 
operation of a system. Examples for expenses are fuel cost (FC) 
and operating and maintenance costs (OMC) [25]. The series of an¬ 
nual costs associated with the carrying charges (CQ) and expenses 
(FCj and OMCf for the jth year of the system operation is not uni¬ 
form. In general, carrying charges decrease while fuel costs in¬ 
crease with increasing years of operation. Therefore, in economic 
analysis of the proposed micro-CHP system, levelized values for 
costs namely CC L , FC L , OMC L and TRR L were used. In this analysis, 
the payback period for additional investment of the micro-CHP 
system (compared to the reference building) denoted as PBP was 
obtained using the levelized total revenue requirement ( TRR L ), lev¬ 
elized cost of energy saving (CES L ), levelized income from electric¬ 
ity sold to the local grid (IES L ) and levelized reduction in the 
associated social cost for carbon dioxide ( CDERC L ) from the follow¬ 
ing expression: 


+ 


/ Qhws 

\ R hws 


Qr + Qch + 



(10) PBP 


TRRl 


■ CDEF 


CES l + IES l + CDERC l 


(13) 
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Procedures for the estimation of TRR L , CES L , IES L and CDERC L are pre¬ 
sented in the following sections. 


3.6.1. Calculation of total levelized revenue requirement, TRR L 

According to the procedure given in [25], the levelized total rev¬ 
enue requirement of an energy system was obtained. In this paper, 
this procedure was employed and adjusted based on feature of the 
proposed micro-CHP system. It should be noted that tax and insur¬ 
ance rates were ignored in this economic analysis as these eco¬ 
nomic parameters were not significant or existed for small 
energy systems. For the levelized total revenue requirement, the 
following can be given [25], 



where CRE is capital recovery factor obtained from the following 
expression [19]: 


CRF 


ieffj 1 + ie//) BL 

(1+Wf -1 



In Eqs. (14) and (15), i eff and BL are effective interest rates (assumed 
as 12%) and booked life of the micro-CHP system (assumed as 
20 years), respectively. In Eq. (14), k is a factor obtained from infla¬ 
tion rate (escalation factor), r (assumed as 6%) and effective interest 
rate (i e //) as follows [20], 


k 


1 + r 

1 + ieff 



TRRj is the total revenue requirement associated to the additional 
investment of the micro-CHP system (compared to investment of 
the reference building) for thejth year of the system operation. TRRj 
for the solar dish Stirling micro-CHP system is simplified and calcu¬ 
lated from the following expression, 

TRRj = TCRj + ROIj + OMCj ( 17 ) 


where TCRj , ROIj and OMCj are total capital recovery, return on 
investment and operating and maintenance costs for jth year of 
the system operation, respectively. These parameters were deter¬ 
mined based on the procedure given in following sections. 


3.6.1.1. Calculating total capital recovery, TCRj. TCRj is obtained 
based on the total depreciable investment (TDJ) and booked life 
of the system (BL = 20) as follows (from [25], with slight 
modification), 

TCRj = ™ for j= 1 , 2 ,..., BL (18) 


3.6.12. Calculating return on Investment, ROIj. Return on investment 
for jth year of the system operation ROIj was obtained based on 
expression given in [25] with slight modification as follows, 

ROIj = BBYj ■ i for j = 1,2,..., BL (19) 

where BYj is the balance at the beginning of the jth year and iis the 
annual rate of return for investment, which was assumed 15% in 
this study. With slight modification from [25], the following can 
be written 

BBYj 

where TNI and TDI are total net investment and total depreciable 
investment, respectively. For simplicity, it was assumed that total 
net investment was depreciated in booked life of the system; there¬ 
fore: TDI = TNI. 


TNI 


for j = 1 




BBYj_ ! -m for j = 2,3,. ..,BL 


( 20 ) 


Calculation procedure for TNI of large energy systems is compli¬ 
cated and is obtained from several investment parameters as pre¬ 
sented in [25]. Fortunately, the present system was a small energy 
system and its investment was simply composed of the purchased 
equipment cost, engineering and construction cost, contingency 
and other costs. Other components of the total net investment such 
as start-up cost, piping and instrument cost, engineering cost, land 
cost and so on were not significant or not exist in the case of the 
micro-CHP system. Therefore, only the additional cost which was 
invested to purchase and employ the solar dish Stirling system 
was considered here. Total purchased equipment cost (TPC) was 
composed of the purchased cost of the solar dish (including collec¬ 
tor and receiver), free piston Stirling engine, heat exchanger#2 and 
domestic hot water heater#l (D.W.H.l was used to recover the re¬ 
jected heat from the cold end of the Stirling engine-see Fig. 1 ). Fur¬ 
thermore, when D.W.H.l was employed in the CHP system, it 
reduced the required heating capacity of the hot water boiler; 
hence, a cheaper hot water boiler should be purchased. Therefore, 
difference between the purchased cost of the hot water boiler of 
the reference building and that of micro-CHP system should be 
considered in calculation of the total purchased cost of CHP equip¬ 
ment (TPC). Thus, 

TPC = Ccolec + Cresiv + Qe + Ov.E2 + ( Cb t REF ~ Cb,CHP ) ( 21 ) 

Ccolec. Cresivi Cse. Ch.e2> Qj,ree and Ct, t cHP are purchase equipment 
costs for solar collector, solar receiver, free piston Stirling engine, 
heat exchanger #2, hot water boiler used in the reference building 
and hot water boiler used in the CHP system, respectively. These 
data were obtained from vendor price list of the equipment. 


colec — 400 Acolec US $ 

(22a) 

RESIV = 0.250(DM) US $ 

(22b) 

se = 500. W S £ US $ 

(22c) 

h.£2 = 6.1 745Qh. £ + 134.5 US $ 

(22d) 

= 10.43Qfa +94.572 US $ 

(22e) 


where Qh.e and Qb are the heat duty of heat exchanger (kW) and 
heating capacity of the hot water boiler (kW), respectively. Contin¬ 
gency, piping and installation and other costs were considered 10%, 
16% and 3% of the total purchased cost of equipments (TPC); there¬ 
fore, for total net investment, the following can be assumed 

TNI = 1.29TPC (23) 


3.6.1.3. Calculating Operating and Maintenance Cost, OMCj. In this 
section, the procedure for evaluating annual cost of operating 
and maintenance for jth year of the system operation is presented. 
Since additional cost of the micro-CHP system was to be compared 
with the reference building using conventional HVAC system, only 
the associated operating and maintenance cost of the dish Stirling 
part was considered in this analysis. Based on the expression given 
in [25], operating and maintenance cost of a system for jth year of 
the system operation could be obtained as, 

OMCj = ( koM cd ~ ,(i 0Mc) \ . 0MCi j or j = 2 ,3 (24) 

VI— I<omc J 

where OMCi is the operating and maintenance cost in the first year 
of the system operation. It is assumed to be 5% of purchased cost for 
the CHP equipment. k 0MC is a factor that is obtained from the similar 
expression to Eq. (16) but with corresponding escalation factor for 
the cost of operating and maintenance, r 0 Mc (assumed as 5%). 
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3.6.2. Calculating levelized cost of energy saving, CES L 

The levelized cost of energy saving was evaluated based on the 
magnitude of primary energy saving and levelized energy price 
(fuel cost) for each energy carrier including natural gas and elec¬ 
tricity. Therefore, 

CES l = FC Lng + FC Le (25) 

where FC Lng and FC Le are levelized annual fuel cost (energy price) of 
natural gas and electricity, respectively. The levelized fuel cost for 
each energy carrier including natural gas and electricity was de¬ 
fined as [25], 

FQ = CRF ■ • FCi (26) 

where k FC is a factor that is obtained from the similar expression to 
Eq. (16) but with corresponding escalation factor for the cost of fuel 
cost, r FC (assumed as 6%). FQ is annual fuel cost in the first year of 
the system operation, evaluated from the following expression, 

EC ! = (c fNG ■ PES ng ) + (c fE ■ PES e ) (27) 

where PES N g and PES e are primary energy saving for natural gas and 
electricity and C/ NG and cy £ are unit costs of natural gas and electric¬ 
ity, respectively. These unit costs were evaluated based on OECD 
countries’ market as 0.0126$ kW h -1 and 0.015$ kW h -1 for natural 
gas and electricity, respectively, as presented in Table 13. 


3.6.3. Calculating levelized income from energy sold to local grid , IES L 
The proposed dish Stirling system had flexibility to buy and sell 
electric energy from/to the local grid. As mentioned previously, 
when the output capacity of the solar dish Stirling engine was low¬ 
er than the total electric demand, difference between the actual de¬ 
mand and power generation was provided by local grid. On the 
contrary, when electricity demand was lower than the output 
power, the excessive electric energy was sold to the local grid. 
Therefore, in this case, there was income from the electricity pur¬ 
chased to the grid. In calculation of the system payback period, the 
levelized value for this income should be evaluated as predicted in 
Eq. (13). For IES L} the following can be given, 


IES l = CRF ■ 


A(i-fcf) 
( (i - h) 


■IESi 



where k E is a factor obtained from the similar expression to Eq. (16) 
but with corresponding escalation factor for the cost of electricity, r E 
(assumed 6%). JFSi is the annual income from electric energy sold to 
the local grid in the first year of the system operation evaluated 
from the following expression, 

IESi = c E ■ EESG (29) 

where c E and EESG are unit costs of electric energy sold to the grid 
(assumed as 0.02US$ kW h -1 ) and annual electric energy sold to the 
grid in kW h, respectively. For EESG, the following can be given, 

^ 8760 h 

EESG = / Esoid dt (30) 

Jo 

where E sold is obtained from Eq. (3b). 


3.6.4. Calculating levelized carbon dioxide emission reduction sost , 
CDERC l 

Carbon dioxide emission imposes the social cost associated 
with health and environmental expenditure. When the solar dish 
Stirling CHP is implemented, due to reduction in greenhouse gas 
emission, there is a kind of monetary saving in the social scale. 
In calculation of the payback period for return of investment of 
the CHP system, such a kind of monetary saving should be consid¬ 


ered as predicted by Eq. (13). Therefore, a levelized value for the 
associated cost of carbon dioxide emission reduction ( CDERC L ) is 
calculated as follows, 


CDERC l 


CRE 


I<CDERC( 1 — l< 


.BL 
CDERC, 


(1 — k CDERC ) 


• CDERCi 



where k CDERC is a factor obtained from the similar expression to Eq. 
(16) but with corresponding escalation factor for the cost of elec¬ 
tricity, r cderc (assumed as 6%). CDERCi is the annual cost associated 
with carbon dioxide emission reduction in the first year of the sys¬ 
tem operation. CDERCi is evaluated from the following expression, 

CDERCi = c e nv. ■ CDER (32) 

where CDER is annual reduction of the carbon dioxide emission as 
evaluated in Section 3.5. c env . is the environmental cost associated 
with each kilogram of the carbon dioxide emission. In this paper, 
it was considered to be 0.02086 $ kg -1 [26-28]. 


4. Methods 

As mentioned previously, three scenarios were considered in 
analyzing the proposed dish Stirling micro-CHP system. First of 
all, the proposed system as shown in Fig. 1 was analyzed in five cit¬ 
ies with diverse climatic conditions based on the first law of ther¬ 
modynamic using the method mentioned in Sections 3.1-3.3. 
Therefore, load estimations, material and energy balance (as in Ta¬ 
ble 10), first law analysis of the solar dish Stirling engine and final¬ 
ly first law analysis of the entire system were performed in this 
step. 

In the second step, the primary energy saving, reduction on 
carbon dioxide emission and payback period for additional invest¬ 
ment of the CHP system were analyzed and compared to the ref¬ 
erence building with conventional HVAC system in the 
aforementioned five locations. In this section, the proposed CHP 
system was analyzed based on the models presented in Sections 
3.4-3.6. Therefore, values of the annual primary energy saving 
(PES), annual carbon dioxide emission reduction (CDER) and pay¬ 
back period (PBP) for each scenario were obtained for each 
location. 

In the third step, PES, CDER and PBP were considered three key 
criteria for selection and sizing of the suitable CHP system in each 
site. In this regard, for each city, a proper scenario for decision 
making on the size of the solar dish Stirling engine was introduced 
by a process of decision making, in which PES, CDER and PBP were 
criteria of decision making. In this case, decision making was per¬ 
formed using the TOPSIS decision making method. The ideal opti¬ 
mal configuration is obtained in a scenario having maximum PES, 
maximum CDER and minimum PBP. It is clear that each scenario 
might have better values for some criteria but worse values for 
other criteria as PES, CDER and PBP were conflicting criteria. The 
ideal scenario had minimum values for payback period ( PBP idea i ) 
and maximum values for the carbon dioxide emission reduction 
and primary energy saving ( CDER idea i and PES idea i) among three sce¬ 
narios. The worst scenario was a scenario with maximum values 
for payback period ( PBP non _ idea i ) and minimum values for the car¬ 
bon dioxide emission reduction and primary energy saving (CDFR- 
non-ideai and PES non _ idea i). In the TOPSIS decision making [29,30], 
deviation of each scenario from the ideal and non-ideal scenarios 
was evaluated and a proper scenario with minimum deviation 
from the ideal scenario and maximum deviation from the non¬ 
ideal solution was introduced. In the following section, basics of 
TOPSIS decision making are introduced. 
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4.1. Decision making 

In a multi-criteria problem, a process of decision making is re¬ 
quired for selecting the desired option from available options. 
There are several methods for decision making process. Since 
dimension of various criteria in a multi-criteria problem might 
be different, prior to any decision making, dimensions and scales 
of criteria should be unified. In this regard, criteria vectors should 
be non-dimensionalized before decision making. There are some 
methods of non-dimensionalization utilized in decision making 
including linear non-dimensionalization, Euclidian non-dimen- 
sionalization and fuzzy non-dimensionalization. As Euclidian 
non-dimensionalization was employed in TOPSIS method [29,30], 
Euclidian non-dimensionalization is discussed as follows, 


d ,r = ~ nnon-Ueaf (34b) 

where m denotes number of criteria. In continuing the TOPSIS 
method, a parameter denoted by c, can be defined as follows, 


In TOPSIS method, a solution with maximum c, is selected as a de¬ 
sired final solution; therefore, if ifi nah is an index for the final se¬ 
lected solution, then, 

ifinai = i e max(Cj)i = 1,2,...,m (36) 

5. Results and discussion 


4.1.1. Euclidian non-dimensionalization 

In this method, a non-dimensionalized criterion, F'\ is defined 
as, 



where F, denotes value of each criteria (Fi = PES, F 2 = CDER and 
F 3 = PBP ) and m is total number of criteria (m = 3 in this study). 


4.2. TOPSIS decision making method 


An ideal solution is the solution with the best value for each cri¬ 
terion and the non-ideal solution has the worst value for each cri¬ 
terion among the values obtained for that criterion in all options 
(scenarios in this study). It is obvious that, in the multi-criteria 
problem with conflicting criteria as in the present study’s case, it 
is impossible for any scenario to have the best values for all the cri¬ 
teria. Therefore, the ideal option (scenario) does not exist in real 
world. In the TOPSIS method, after Euclidian non-dimensionaliza¬ 
tion of all the criteria, an ideal criterion and a non-ideal criterion 
are specified through comparison of values of each criteria (F,) 
for all scenarios. Therefore, the best values of all the criteria are 
dedicated to the ideal solution (F, >IC / e a/) and the worst values are 
dedicated to the non-ideal solution ((F i non _ lde a/))- Then, all the crite¬ 
ria are non-dimensionalized based on Eq. (30) and special dis¬ 
tances from the ideal and non-ideal solutions denoted by dt and 
df in the non-dimensionalized criteria space are determined as 
follows, 



Yl7 F i- F im,y 


(34a) 


5.1. Monthly variation of key parameters 

Figs. 7-10 illustrate monthly variations of the Stirling engine’s 
electric efficiency, overall efficiency of the solar dish Stirling CHP 
system, primary energy saving (PES) and carbon dioxide emissions 
reduction (CDER) for Scenario I (sizing of the primary mover based 
on the minimum annual electricity demand) at various locations, 
respectively. Fig. 7 demonstrates monthly variation of thermal effi¬ 
ciency of the solar Stirling engine (Scenario I). It indicates that 
maximum thermal efficiency of the solar Stirling engine was 
34.5% obtained in Tabriz in May. On the other hand, minimum effi¬ 
ciency was 25.9% in Rasht in December. These values could be 
attributed to the variation of direct normal solar irradiation which 
was at the highest level in summer and the lowest level in winter. 
Moreover, low ambient temperature in Tabriz caused higher Car¬ 
not efficiency and subsequently more real electric efficiency. 

Fig. 8 represents monthly variation of the overall energy effi¬ 
ciency of the micro-CHP system for the system sized based on Sce¬ 
nario I. The highest and lowest efficiency rates were 78.4% and 
53.1%, which were associated with Tabriz in March and Bandar Ab¬ 
bas in August, respectively. Considering Eqs. (16) and (18) pre¬ 
sented in Table 10, since the maximum heat recovery (Qr) 
occurred in cold months, it is anticipated to have the maximum 
efficiency in these months. On account of having the lowest hot 
water demand and consequently the lowest heat recovery in Ban¬ 
dar Abbas, it can be predicted that the minimum overall efficiency 
can be obtained in this city. 

Fig. 9 demonstrates percentage of primary energy saving (PES) 
for the solar dish Stirling CHP system in comparison to the refer¬ 
ence building that did not utilize the CHP system. The PES percent¬ 
age highly depends on four factors: solar direct normal irradiation, 


Tehran Tabriz — — Bandarabbas — • -Rasht .Yazd 



Fig. 7. Monthly Stirling engine efficiency for the selected cities. 
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Tehran-Tabriz — — Bandarabbas — • -Rasht .Yazd 



Tehran Tabriz — — Bandarabbas — • -Rasht .Yazd 



Tehran Tabriz — — Bandarabbas — • - Rasht.Yazd 



Month 

Fig. 10. Monthly percentage of carbon dioxide emissions reduction for the selected cities. 


ambient temperature, number of hours with solar radiation per 
day and excess electricity produced by the Stirling engine. How¬ 
ever, the most prevalent factor is the solar direct normal irradia¬ 
tion. As shown in Fig. 9, the maximum primary energy saving 
(16.53% in April) was acquired in the case that this CHP system 
was utilized in Rasht since the excessive electricity produced by 
the Stirling engine (E ex ) in this location was the highest among 
other locations, which led to higher primary energy saving. It is 
essential to mention that the minimum primary energy saving 


was associated to Tabriz in January with value of 3.25% which 
was influenced by four cited factors. 

Fig. 10 illustrates monthly variation of the percentage of carbon 
dioxide emission reduction which was gained when the solar dish 
Stirling system CHP (Scenario I) was utilized. Similar to the trend of 
the PES variation, maximum and minimum percentage of carbon 
dioxide emissions reduction belonged to Rasht in April and Tabriz 
in January, respectively (22.8% and 3.9%). It can be concluded that 
the variation trend of carbon dioxide emission reduction was 
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similar to that of the primary energy saving because emission of 
carbon dioxide was almost influenced by consumption of electric¬ 
ity and natural gas in the site, amount of heat recovery from the 
Stirling engine and amount of electric energy sold to the grid. 

5.2. Annual values of key parameters 

Figs. 11-15 demonstrate annual values of the average thermal 
efficiency of the Stirling engine, average overall efficiency of the so¬ 
lar dish Stirling CHP system, primary energy saving (PES), carbon 
dioxide emissions reduction and values of the payback period pre¬ 
sented in three scenarios for the selected locations, respectively. 
These data were presented in order to have better insights regard¬ 
ing overall operation of the proposed CHP system. 

Fig. 11 shows average annual thermal efficiency of the Stirling 
engine efficiency obtained at different locations. Since thermal effi¬ 
ciency of the Stirling engine in this analysis did not depend on the 
size of engine, value of this efficiency at each site was the same in 
the three scenarios. Therefore, Fig. 11 does not separately show 
Scenarios I, II and III. The highest value of average annual efficiency 
of the Stirling engine efficiency belonged to Bandar Abbas (32.7%) 
and the lowest value was associated to Rasht (31.5%). On the ac¬ 
count of considering the real engine efficiency as the percentage 
of Carnot efficiency in this study, it is anticipated that the most 
important factor affecting thermal efficiency could be difference 


■ scenario I ■ scenario II scenario III ° 

ro 



Tehran Tabriz Bandar-Abbas Rasht Yazd 


City 

Fig. 13. Annual primary energy saving for the selected cities. 

between the temperature of heat source and heat sink in Stirling 
engine operation. Therefore, due to having the maximum differ¬ 
ence between these two temperatures in Bandar Abbas, it can be 
concluded that this city had the highest Stirling engine efficiency. 

Average annual overall efficiency of the solar dish Stirling CHP 
system for three scenarios and five locations is presented in 
Fig. 12. This parameter depended on the series of factors such as 
energy demands of the proposed building, solar irradiation, output 
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Fig. 11. Annual values of Stirling engine efficiency for the selected cities. 
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Fig. 12. Annual values of the overall efficiency of the CHP-Stirling system for the selected cities. 
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Fig. 14. Annual carbon dioxide emissions reduction for the selected cities. 
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Fig. 15. Values of the payback period of different scenarios obtained for the 
selected cities. 


electric power of the Stirling engine and amount of heat recovery 
at the cold-end of the Stirling engine. The figure shows that Tabriz 
had the maximum overall efficiency (71.7%, 71.96% and 74.08% for 
Scenarios I, II and III, respectively) among the selected locations. 
This could be due to the fact that, in Tabriz, the inlet tap water tem¬ 
perature to the D.W.H#1 was in the lowest value; therefore, this 
city had superior capability of absorbing the heat recovered from 
the cold-end of the solar Stirling engine. On the other side, due 
to a similar reason in Bandar Abbas with the highest inlet tap water 
to the D.W.H #1, the lowest average annual overall efficiency oc¬ 
curred (54.87%, 55.23% and 55.44% for Scenarios I, II and III, respec¬ 
tively). It should be cited that Tehran and Yazd had very little 
difference in the amount of system’s overall efficiency (67.5% and 
67.2%, respectively) because of being similar in the building energy 
demands, solar irradiation and amount of heat recovered in both 
cities. 

Furthermore, Fig. 12 indicates that, in all the locations, the CHP 
system that was sized based on Scenario III had the highest value 
for the average overall thermal efficiency of the CF1P system, which 
was for the sake of the fact that, in this scenario, size of Stirling sys¬ 
tem was higher than the corresponding size of the Stirling engine 
in Scenarios I and II. Therefore, by considering reduction of param¬ 
eters in denominator of Eq. (4), the overall efficiency in Scenario III 
was higher than the corresponding overall efficiency in Scenarios I 
and II. 

Fig. 13 illustrates total amount of annual primary energy saving 
for three scenarios and five locations. The primary energy saving 
for the presented system hinged on some factors such as amount 
of solar irradiation, annual number of hours with solar irradiation, 
amount of heat recovered from the cold-end of the Stirling engine 
and excessive electric power generated by the Stirling engine. The 
presented figure shows that the highest annual primary energy 


saving belonged to Bandar Abbas. It revealed that, in Bandar Abbas, 
higher solar irradiation condition was more dominant and influen¬ 
tial factor than the amount of heat recovered at the cold-end of the 
Stirling engine. The minimum annual primary energy saving was 
related to Rasht which had the worst annual solar irradiation con¬ 
dition among the proposed locations. Moreover, the mild outdoor 
winter design temperature in Rasht led to decreasing the amount 
of heat recovery in comparison to the cold cities. This is another 
reason for causing lower primary energy saving in Rasht in com¬ 
parison to other locations. Furthermore, as is evident from 
Fig. 13, the highest primary energy saving was in Scenario III for 
all the locations, which was on the ground that, in Scenario III, 
there was a Stirling engine with greater capacity in comparison 
to other two scenarios for all the locations. Therefore, in Scenario 
III, there was lower purchased electric energy from the grid and 
lower natural gas consumption (because of higher heat recovery 
at the cold-end of engine); hence, PES was higher in this scenario. 

Fig. 14 demonstrates annual amount of carbon dioxide emission 
reduction for three scenarios and five locations. As carbon dioxide 
emission was directly proportional to the amount of primary en¬ 
ergy consumption, due to a similar reason to the previously ex¬ 
pressed one for the PES, in terms of annual carbon dioxide 
emissions reduction among the selected cities, Rasht and Yazd 
had a lowest and greatest promising potentials in reducing the 
amount of carbon dioxide emission in Scenarios II and III respec¬ 
tively. In Scenario I, carbon dioxide emission reduction potential 
of Tabriz was a little lower than that of Rasht, due to having unfa¬ 
vorable solar irradiation in Tabriz as a result of the cloudy sky in 
most times of the year. In a similar trend to the primary energy 
saving, value of carbon dioxide emission reduction for Scenario 
III was greater than the corresponding values of Scenarios I and II 
for all the locations. 

Figs. 13 and 14 reveal that, in terms of the primary energy sav¬ 
ing and carbon dioxide emission reduction, Tabriz and Rasht were 
two similar worst cities which could be attributed to the similarity 
of direct normal solar irradiation, ambient temperature, electricity 
demand of the building and amount of power produced by the Stir¬ 
ling engine. 

Payback period for additional investment of the solar dish Stir¬ 
ling CHP system as a very important factor for assessment of the 
proposed CHP system for three scenarios and five locations is pre¬ 
sented in Fig. 15. According to the data presented in Fig. 15, it can 
be seen that Yazd had the best economic condition for utilization of 
the solar dish Stirling CHP system. Values of payback period for 
this city in Scenarios I, II and III were 6.2, 8.8 and 10.5 years, 
respectively. With respect to Eq. (13), the payback period de¬ 
pended on four factors including levelized values for the total rev¬ 
enue requirement (TRR L ), cost of energy saving ( CES L ), income from 
electricity sold to the local grid (IES L ) and carbon dioxide emission 
reduction cost ( CDERC L ). TRR L was entirely affected by the cost of 
system (included investment and O&M costs) with respect to the 
reference building. Other three mentioned parameters were re¬ 
lated to energy and environmental operation of the system. There¬ 
fore, in the case of having good solar irradiation in Yazd, it is 
expected to have one of the most primary energy saving costs 
and carbon dioxide emission reduction costs, which consequently 
causes the lowest payback period for this city. In the case of Rasht 
due to having the worst solar irradiation, the values of payback 
period for all three scenarios were much higher than other cities. 
It should be mentioned that, because of having undesirable solar 
irradiation in Rasht, higher size of dish Stirling system was re¬ 
quired in order to cut down the value of payback period. This fact 
is clear as shown in Fig. 15 for Scenario III in Rasht. On the other 
side, sizing of the prime mover based on maximum annual electric¬ 
ity demand (Scenario II) in cities which had fairly good solar irra¬ 
diation (Tehran, Tabriz and Yazd) led to increase in the value of 
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Table 14 

Decision making based on the TOPSIS method for five locations. 


Scenario 

PES (MW h year” 1 ) 



CDER (Ton year 




Payback period (year) 



d Index (TOPSIS) 




Tehran 

Tabriz 

Bandar- 

abbas 

Rasht 

Yazd 

Tehran 

Tabriz 

Bandar- 

abbas 

Rasht 

Yazd 

Tehran 

Tabriz 

Bandar 

abbas 

Rasht 

Yazd 

Tehran 

Tabriz 

Bandar 

abbas 

Rasht 

Yazd 

I 

71.5 

62.1 

116.7 

55.3 

77.8 

9.9 

8.5 

18.2 

8.7 

10.7 

6.6 

7.8 

13.0 

27.8 

6.2 

0.622 

0.621 

0.425 

0.388 

0.623 

II 

87.9 

72.2 

111.6 

41.9 

93.4 

13.0 

10.4 

17.3 

6.3 

13.6 

9.5 

10.4 

12.3 

20.3 

8.8 

0.499 

0.529 

0.453 

0.406 

0.510 

III 

145.2 

137.4 

133.8 

72.9 

165.3 

23.2 

22.0 

21.2 

11.8 

26.5 

12.1 

12.5 

11.2 

13.0 

10.5 

0.227 

0.219 

0.397 

0.546 

0.225 

Ideal 

145.2 

137.4 

133.8 

72.9 

165.3 

23.2 

22.0 

21.2 

11.8 

26.5 

6.6 

7.8 

11.2 

13.0 

6.2 

1.0 

1.0 

1.0 

1.0 

1.0 

Non¬ 

ideal 

71.5 

62.1 

111.6 

41.9 

77.8 

9.9 

8.5 

17.3 

6.3 

10.7 

12.1 

12.5 

13.0 

27.8 

10.5 

0.0 

0.0 

0.0 

0.0 

0.0 


Table 15 

Final specifications of the selected solar dish Stirling micro- 

■CHP systems at five locations. 







Parameter 

Tehran 


Tabriz 


Bandar Abbas 

Rasht 


Yazd 



Actual 

value 

Value per 
m 2 of 
building 

Actual 

value 

Value per 
m 2 of 
building 

Actual 

value 

Value per 
m 2 of 
building 

Actual 

value 

Value per 
m 2 of 
building 

Actual 

value 

Value per 
m 2 of 
building 

Output power (kW) 

4.00 

0.0040 

3.58 

0.0036 

9.60 

0.0096 

8.00 

0.0080 

3.80 

0.0038 

Dish diameter (m) 

5.05 

0.0051 

4.72 

0.0047 

7.5 

0.0075 

10.00 

0.0010 

4.88 

0.0049 

Annul average thermal 

31.81 

N.A 

31.68 

N.A 

32.75 

N.A 

31.52 

N.A 

32.47 

N.A 

efficiency of the Stirling 
engine (%) 

Overall annual average 

67.56 

N.A 

71.72 

N.A 

55.23 

N.A 

67.44 

N.A 

67.26 

N.A 

efficiency of the CHP 
system (%) 

Annual primary energy 

71.475 

0.0715 

62.125 

0.0621 

111.634 

0.1116 

72.928 

0.0729 

77.792 

0.0778 

saving (MW h year -1 ) 
Annual carbon dioxide 

9.9 

0.0099 

8.5 

0.0085 

17.3 

0.0173 

11.8 

0.0118 

10.7 

0.0107 

emission reduction 
(Ton year -1) 

Total net investment (US $) 

25510.0 

25.510 

22415.2 

22.415 

54569.9 

54.570 

103139 

103.139 

23950.0 

23.950 

Levelized total revenue 

1845.4 

1.845 

1621.5 

1.622 

3947.6 

3.948 

7461.2 

7.461 

1732.6 

1.733 

requirement ( TRR L ) (US $) 
Levelized cost of energy 

124.5 

0.124 

107.5 

0.108 

204.3 

0.204 

136.6 

0.137 

135.1 

0.135 

saving ( CES L ) (US $) 
Levelized income from 

2.1 

0.002 

0.0 

0.000 

70.8 

0.071 

304.7 

0.305 

1.2 

0.001 

electricity sold to the grid 
(. IES l ) (US $) 

Levelized cost of carbon 

26.2 

0.026 

22.4 

0.022 

45.6 

0.046 

31.2 

0.031 

28.3 

0.028 

dioxide emission 
reduction ( CRERC L ) (US $) 
Payback period (year) 

6.6 

N.A 

7.8 

N.A 

12.3 

N.A 

13.0 

N.A 

6.2 

N.A 


payback period. In Tehran, Tabriz and Yazd, Scenario I was more 
desirable and Scenario III was the worst scenario from the view¬ 
point of the payback period. In Rasht and Bandar Abbas, the best 
scenario was Scenarios III and I was the worst one considering 
the payback period of the system investment. One common feature 
of these two cities is humid weather. Therefore, it can be said that, 
for cities with humid weather, Scenario III is a better option than 
the other two. 

5.3. Specifying the best scenario for sizing solar dish Stirling micro-CHP 
system 

In previous Sections 5.1 and 5.2, basic features of the proposed 
solar dish Stirling micro-CHP system were assessed based on key 
parameters of the system including thermal efficiencies, primary 
energy saving (PES), carbon dioxide emission reduction (CDER) 
and payback period (PBP) for additional investment of the CHP sys¬ 
tem. In this section, the best scenario for each location is suggested 
based on three criteria of PES, CDER and PBP. In this regard, ideal 
and non-ideal scenarios are defined based on the best and worst 
values of three criteria, as described in Section 4. Therefore, for 
each location, the best scenario was identified using the TOPSIS 


decision making method, as described in Section 4.2. Table 14 rep¬ 
resents values of three criteria for Scenarios I, II and III, both ideal 
and non-ideal, in five locations. Furthermore, this table shows q in¬ 
dex (based on Eq. (35) as used in TOPSIS decision making method) 
for each scenario and location. As in the TOPSIS method, an option 
with the highest value of q was desired (according to Eq. (36)); for 
each location, a scenario which had the highest q was selected as 
the desired scenario for sizing the micro-CHP system. The highest 
q index in Table 14 is specified by bold-underlined fonts. According 
to Table 14, the final recommended scenario for Tehran, Tabriz and 
Yazd was the system sized based on Scenario I (sized based on 
minimal annual electric power demand). For Bandar Abbas, Sce¬ 
nario II in which the system was sized based on the maximum an¬ 
nual value for electric power demand was recommended. For 
Rasht, a constant size of the solar Stirling system based on Scenario 
III was suggested by the TOPSIS decision maker. 

Table 15 summarizes the finally suggested features of the solar 
dish Stirling micro-CHP system for each location. Furthermore, as 
value of the finally recommended system was presented for a spe¬ 
cific building with 960 m 2 residential area in order to suggest some 
thumb rule values for other buildings, values of parameter per 
square meter of the building are also presented in Table 15. 
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6. Conclusion 

This paper presented thermodynamic analysis of a CHP-solar 
dish Stirling engine which was integrated with a conventional 
HVAC system in order to provide a part of heating, hot water de¬ 
mand and electrical loads of a building including electric demand 
of HVAC equipment. Likewise, considerable reduction or even 
omission of the pollution emissions in power generation section 
of the presented system was a promising objective that could de¬ 
crease the emissions produced by burning fossil fuels in power 
plants. Performance of the system was analyzed in diverse climatic 
conditions. In this regard, five cities in Iran including Tehran (with 
temperate and dry weather), Tabriz (with cold and dry weather), 
Bandar Abbas (with hot and humid weather), Rasht (with temper¬ 
ate and humid weather) and Yazd (with hot and dry weather) were 
considered. Among these cities, from the view point of overall en¬ 
ergy efficiency, Tabriz had the highest overall efficiency but, in 
terms of annual primary energy saving, annual carbon dioxide 
emissions reduction and annual Stirling engine efficiency, Bandar 
Abbas had better performance. It was shown that this system 
had a great potential for reducing primary energy saving and car¬ 
bon dioxide emission reduction. Economic analysis of the system 
indicated that, in Tehran, Tabriz and Yazd, the payback period for 
the investment of the CHP system was almost less than half of 
the system life cycle (in Scenario I, it was almost one third of the 
system life time). In Bandar Abbas, the payback period was a little 
greater than the half-life of the system. However, in Rasht, due to 
the great number of times with cloudy weather, the payback per¬ 
iod was not acceptable for Scenarios I and II and the payback per¬ 
iod was a little greater than the half-life of the system for Scenario 
III. Since economic criterion showed conflicting behaviors with the 
primary energy saving and carbon dioxide reduction criteria, the 
best system configuration was specified through a decision making 
tool called the TOPSIS decision-maker. It was shown that, in cities 
with dry weather such as Tehran, Tabriz and Yazd, the best sce¬ 
nario in sizing the solar dish Stirling system was Scenario I, in 
which the system was sized based on the minimum annual electric 
power demand. In Bandar Abbas with hot and humid weather, the 
best scenario was the scenario in which the system was sized 
based on the maximum annual electric power requirement (Sce¬ 
nario II). In Rasht with temperate and humid weather, the fixed 
size of the Stirling engine (Scenario III) was the best choice. This 
analysis showed a useful methodology and tool in energetic, envi¬ 
ronmental and economic assessment of each residential CHP or 
CCHP system and could be applied for assessment and size evalu¬ 
ation of the prime mover for the CHP or CCHP systems in each cli¬ 
matic and solar irradiation condition. 
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